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Abstract 

The one- and the two-particle propagators for an infinite non-interacting Fermi system are stud- 
ied as functions of space-time coordinates. Their behaviour at the origin and in the asymptotic 
region is discussed, as is their scaling in the Fermi momentum. Both propagators are shown to have 
a divergence at equal times. The impact of the interaction among the fermions on their momentum 
distribution, on their pair correlation function and, hence, on the Coulomb sum rule is explored 
using a phenomenological model. Finally the problem of how the confinement is reflected in the 
momentum distribution of the system's constituents is briefly addressed. 

PACS numbers: 24.10.-i, 24.10.Cn, 25.30.-c 
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I. INTRODUCTION 



In this study we derive expressions in space-time for the one- and two-fermion propaga- 
tors, considering the simple case of a non-relativistic, non-interacting infinite Fermi system 
where these quantities in energy-momentum space are well-known. To obtain the corre- 
sponding space-time results is not an entirely trivial job and we were unable to find any 
detailed discussion of these quantities in the literature. The motivation for the present 
study comes principally from the desire to investigate the roles played by correlations, in 
particular short-range correlations, in nuclear matter and finite nuclei, although our study 
applies as well to other fields of many-body physics. From this perspective the present work 
should be viewed as a first step in the direction of treating more complex systems: we present 
arguments later that the correlations among the constituents of a Fermi |lj system are best 
appreciated in space-time, especially so when compared with the non-dynamically-correlated 
situation. We shall illustrate this point using a model for the dynamical correlations which 
modifies the non-interacting step-function momentum distribution at the Fermi surface. 

Two further items are addressed because of their significance for a non-interacting system. 
The first relates to the propagators' second kind scaling property (independence of the Fermi 
momentum) ,|5|: we prove that, providing an appropriate rescaling of space and time 

is performed, both the one-fermion and two-fermion propagators do scale in kp, the latter 
except at equal times. Indeed one finds that when the quantum field theory (QFT) is applied 
to a finite-density many-body system a divergence occurs if the propagators are evaluated at 
equal times. One also finds that both the one- and two-particle propagators become purely 
imaginary at equal times. 

The second item concerns the problem of the impact of confinement on the non-interacting 
fermion momentum distribution (and hence on the propagators). 

The present paper is organized as follows. In Sec. [Til we deal with the one-body propa- 
gator G°(x,x'). We observe that, while the hole propagator is well-defined, at equal times 
the particle propagator is not. We derive an analytic expression for the space-time hole 
propagator and discuss its asymptotic behaviour. We find that it vanishes as a power law 
at large \x — x'\ and t — t' as a consequence of the cut the function G°(k, uj) displays in 
the complex o;-plane just above the real axis for < uj < uip, where ujp = k F /2m, with 
kp the Fermi momentum and m the fermion (nucleon) mass (Pailey- Wiener theorem {(]]). 
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We also prove that G°(x, x') divided by the density scales in kp. In Sec. [TTT] we proceed 
to study the two-particle propagator, specifically the density-density correlation taction, 
starting by focusing on the Coulomb Sum Rule (CSR) [7]. To illustrate how such topics can 
be addressed in space-time, we derive the well-known expression for the CSR using the pair 
correlation function (which arises in the present case simply from the Pauli principle) as 
input and performing the integration in the complex coordinate space. In Sec. HVI we obtain 
the space-time expression for the density- density correlation function, usually referred to as 
II (x, x'). We do this partly in terms of the error function and partly in terms of a function 
for which we keep the integral representation, although it may be expressed in terms of the 
Meijer G- functions |8|]. We analyze the asymptotic behaviour of U°(x,x'), compare it with 
that of G°(x, x') and discuss its scaling behaviour in kp. Next we show that the imaginary 
part of the density-density correlation function (a branch of the two-particle propagator), 
unlike the CSR, cannot be expressed directly through the Pauli pair correlation function, 
its only ingredient being the momentum distribution. Furthermore we show that the QFT 
shortcoming previously found in the case of G°(x,x') also affects U°(x,x'). In Sec. [V] we 
introduce a phenomenological momentum distribution, which we employ for computing the 
CSR and discuss the significance of the difference between the latter and the CSR of the 
free Fermi gas. In the Conclusions (Sec. IVIj) we briefly address the problem of how the 
confinement of our system affects the momentum distribution of its constituents, summarize 
our findings and outline a few further important issues we intend to address in future work. 

II. THE ONE-BODY PROPAGATOR 

In this section we deal with the four-dimensional Fourier transform of the well-known 
one-body fermion propagator G°(k,u) in an infinite, homogenoeus, non-interacting system. 
That is, we compute 

G°(x,x>) = [ A e *(«) r *(*-**) + 9 ^~ k ) } , (!) 

where here for simplicity spin indices are suppressed. The frequency integration in Eq. (CQ) 
is easily performed in the complex cu-plane and one can recognize forward (particle) and 
backward (hole) propagation. In Eq. (ffj) the angular integrations are also immediate, and 
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one is left with the expression 

■■kp 



G°(x, x') = -^-( 9(t' - t) [ F dkke-^-V sin(fcr) - 0(t - t') [ dkke"^^ sin(Ax) \ 
2tt r { J J kp J 

where r = \x — x' | and Uk = k 2 / (2m). The integral 



(2) 



»/7T VP 2(t-t') 

dkke-^-^ sin(A;r) = ^ (3) 



:^-*')) ; 

diverges when computed at equal times t = t' for any value of r, and thus at equal times 
the particle propagator in the present framework is ill-defined and cannot be computed 
unless regularized. For instance, this might be achieved in the context of the Wightman 
formulation of QFT (see j^]) which replaces the field function with a distribution. 

We do not dwell here on the problem of the regularization of the particle propagator and 
instead limit our attention to the computation of the hole propagator. This can be done 
analytically and yields (a and (3 are the spin indices which we temporarily reintroduce) 



k F r l-i -,[(^) 2 + a,] 



2 

(k F r) 2 + 2iA t 



. , cos(A;irr) sin(fcpr) 
3i(k F r)+ r + 



k F r 2iA t 



2 (2A t ) 3 / 2 2 

x 



where n = kp /3n 2 is the system's density for spin 1/2 particles (the case we are considering). 
Furthermore, 

is the time difference expressed in inverse Fermi frequency units (the natural choice in the 
present context), ji(kpr) is the spherical Bessel function of order one and the standard 
definition 

erf(z) = r e~ t2 dt (6) 
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for the error function is employed. Note that, but for the overall factor n , the Green function 
G°(x,x') scales in the system's density, i.e., it loses all explicit kp dependence which then 
enters only in defining the scale of the space at any value of A t . 

/,From the asymptotic expansion of the erf(z) (see Appendix |A]) it follows that the 
equal-time limit of G°(x, x') reads 

G Q {xt,x't + ) = i r ^^- r J l {k F r) (7) 



which, when r — > and the spin trace is taken, reduces to the system's density, as it 
should. Notably Eq. (jTJ) also holds for finite A f and large r, as is shown in Appendix [A] 
where it is proven that under these conditions the contribution stemming from the two error 
functions exactly cancels the one arising from the second and third terms inside the square 
brackets in Eq. (j3J). Thus, for fixed A t the propagator in Eq. (J3J) asymptotically behaves 
as cos(kpr) / (kpr) 2 , reflecting the square root singularity k = \f2muj — it of G°(k,u) in 
the complex fc-plane. Indeed as illustrated in [6], which deals with the theory of potential 
scattering, the structure of the Fourier transform basically implies that the transform of a 
singularity is an asymptotic behaviour. Analogously G (k,u) displays a [At] -1 behaviour 
at large times for fixed r, since in the complex cu-plane the propagator has a simple pole for 
fixed k. 

The above findings are borne out by the results shown in Figs. [H [2] and [3X where one 
sees the behaviour of the modulus of G°(x, x') versus kpr for a few values of A t . Also shown 
in Figs. [2]and [3] are the real and imaginary parts of G°(x, x') for two finite values of A*. As 
is well-known in (k, to) space, these are connected through a dispersion relation. 

Importantly for Af = the hole propagator becomes purely imaginary and displays an 
oscillatory behaviour versus kpr. This accounts for the zeros in its modulus shown in Fig. [TJ 
From the point of view of the Heisenberg principle A t = corresponds to the maximum 
energy (in fact momentum) uncertainty for the propagating particle. Thus the wave function 
of the latter corresponds to a large superposition of plane waves yielding the striking zeros 
seen in Fig. [TJ For A t small, but not vanishing (Fig. [2]), the propagator starts to develop 
a real part. This also oscillates, but its zeros differ from those of the imaginary part. 
Accordingly, now the particle can be found everywhere in space, but of course with a small 
probability when in the vicinity of the zeros appearing in Fig. [TJ As A t grows (Fig. [3]) so 
does the real part of G°(x, x') and the behaviour of the modulus of the propagator becomes 
smoother and smoother, until for At very large it becomes constant: indeed now, in accord 
with the Heisenberg principle, the wave function of the particle becomes a plane wave. 

III. THE COULOMB SUM RULE 

To pave the way to the more general treatment of the two-particle propagator (of which 
the density-density correlation function is a particular branch) to be discussed in the next 
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section, we once more derive the well-known Coulomb sum rule, although this time through 
a somewhat different technique, namely an integration in coordinate space. For this purpose 
we recall (see [10]) that, using standard quantum mechanics, the CSR can be cast in the 
form 

S(q) =< ^{-m-Qil^o > , (8) 

where > is the system's ground state and p(—q) is the density deviation operator defined 
as 

~p(-q)=p(-q)-<^o\p(-mo> ■ (9) 
It is then a straightforward matter to obtain from Eq. the formula 

s(q) = z + j dSdffe-w*-® < ^U^l(y)^p(y)^a(x)\^ > -f$V6(q) (10) 

which expresses the CSR in terms of the fermion fields in the Schroedinger picture (spin 
indices have been reintroduced) always sticking to the model of a homogeneous system 
enclosed in a large volume V having Z = A/2 charges, with A being the total number of 
spin 1/2 particles. In the second term on the right-hand side of Eq. ffTOj) one recognizes 
the equal-time two-particle propagator. This can be expressed in terms of the correlations 
among the particles in the system and indeed, in our simple case, it reads 

2- 



< * O |*],(^*to*j8(0)*a(£)|tfo > = 



Ur, 



2 V k F r ) 



-g 2 (k F r) (11) 



with r = \x — y\. The function ^g 2 (kpr) is usually referred to as the Pauli pair correlation 
function. In the above the direct and exchange contributions clearly appear. Moreover it 
should be kept in mind that basic in deriving Eq. ( iTTjl have been the ^-functions entering in 
G°(k,u). Our aim here is to show that even small variations of Eq. fTTTT) at short distances 
(and hence of the ^-functions in momentum space) can produce quite dramatic changes in 
the CSR (see Sec.|VD. 

We proceed by inserting Eq. (TTT]) into Eq. (TTOT) . getting for the latter 



S(q) = Z + n 2 V5(q) - n 2 V / dfe 



-tq-r ' 



-n 2 VS(q) 
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(12) 



where the elementary angular integrations have been performed and the term arising from 
the direct piece of Eq. ffTTl) is seen to drop out. 



The integral in Eq. (fl2|) can be computed in the complex z-plane using standard tech- 
niques (see Appendix [B] for details), yielding for the Coulomb sum rule the familiar non- 
relativistic expression 10J 

\z[\4-- h(S-Y) Hq<2k F 
S(q) = \ ^ 4kF 2 V 2fc ^ J (13) 

In connection with the above derivation which is valid for a perfect Fermi gas, it should be 
pointed out that the result in Eq. (1131) stems from the exact cancellation of two contributions, 
as shown in Appendix B. We shall then prove in Section |V] that, as anticipated above, even 
a minor modification induced by interactions among the system's constituents of the 9- 
functions entering in G (k,u) is sufficient to disrupt the cancellation in Eq. (1B9[) : hence 
such a modification induces a sizable change of the CSR for large q. Since, as will be 
shown in Sec. |Vl modifying the ^-function around the value k = k F actually corresponds 
to modifying the pair distribution in Eq. (fTTj) for small distances, this outcome is what one 
should expect and it offers a nice example of how the Fourier transform works. 



IV. THE DENSITY-DENSITY CORRELATION FUNCTION 

In this section we analytically compute and explore in space-time the branch of the 
two-particle propagator usually referred to as the density- density correlation function or 
polarization propagator 11(2, y). The expression for the latter is again well-known in energy- 
momentum space for a non-interacting system, as is the fact that its imaginary part provides 
the inelastic scattering cross section for many types of probes of a many-body system. 

The definition of H(x, 0) is the following 

IT(x,0) =-i< V \T(p H (x)p H (0))\y > , (14) 

where the density deviation operators are in the Heisenberg picture, unlike the case of 
the CSR where they were taken in the Schroedinger picture, i.e., where no T-product was 
introduced. To explore this point further we recast Eq. ([141) in the form 

U(x,0) = -z { <*o|T(*t(»)*«(»)*J(0)*/9(0))l*o> 

- < * \fr a (x)* a {x)\* >< ^o|^(0)^(0)|^ >} • (15) 
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A comparison with Eq. ({TO]) then shows that what enters in the CSR is just the above 
expression, however with the T-product replaced by its equal-time specification and it is 
precisely this quantity which is directly connected with the pair distribution function. As 
we have seen (and as we will discuss in more detail later) the latter crucially affects the 
CSR, namely, the frequency integral of the response. 

The response, however, is not directly expressible in terms of the pair distribution cor- 
relation function, but rather the momentum distribution of the system's constituents (in 
our case a ^-function) enters into its definition. It thus appears that the ^-function should 
be viewed as the fundamental ingredient of both the system's response and of the CSR. 
Indeed in Sec. [V] we shall illustrate how the pair distribution function (and hence the CSR) 
is determined by the momentum distribution. 

To proceed further, observe that with a straightforward application of Wick's theorem to 
Eq. (I15p one obtains for the non-interacting Fermi system 

n°(x, 0) = -2iG°{x, 0)G°{x, 0) . (16) 

To compute this we start from the Fourier transform of its well-known expression in energy- 
momentum space 

7-1; 



n°(ff) = (~ 2i )j ^»(fc)G»(Hg) 

/dk 
—6(\k + q\ - k F )6(k F - k) 



q + u % - + irj q - u % + uj^ - irj _ 

(17) 



where the integration on k° has been performed in the complex plane. Next we take the 
inverse Fourier transform of Eq. ([Tj 



n°(x,0) = J ^ye^n°(q) (18) 

and carry out the go _m tegration in the complex plane: this again distinguishes between 
forward and backward time propagation, corresponding to the two terms generated by the 
T-product. Both of these diverge for t x = in accord with what was previously found for 
the one-fermion propagator. Here we recall that only the particle piece of the latter was 
seen to diverge at equal times. We shall return on this point later on. 
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Furthermore, the second term of U°(x, 0), associated with q < 0, describes, the system's 
response in the time-like domain. One obtains 

n°(z,0) = -2i J J^e«*f j^- 3 6(k F - k)6(\k + q] - k F ) 

x f e - i K+sf^ S (g +c iB+^)«.j(_i | )J . (19) 



After some algebra (see Appendix [C] for details) the above expression can be reduced to a 
single integral (we focus on the first piece, since the second one is immediately derived when 
the first is known): 

U°(x,0)6(t x ) = / dse~ iA * s sm(sk F r) 3l (2A t s) 

+ J dse~ iAts2 sin(sA; F r)i 1 (2A t s(s - l))(s - l) 2 

+ 2A~ / dS sin ( skF ^ e " lAtS2 C0S ( 2A * S ( S - !)) - e ~ 2iAts2 Q cos ( A * s ( s - 2)) 

+ 2~h sin ( AAs - 2) )) ] + ak I ds sm ^^) e_24AiS si ^ 2A ^ 2 ) + *) 

+ cos ^: ) ~ 1 ) , (20) 



4A t A; jF r 

where r = |x| and s = q/kp. As previously anticipated, in the limit of vanishing A t the 
above becomes purely imaginary and diverges. 

To show that this divergence is as severe as the one encountered for the one-particle 
propagator (see Eq. (jHJ)) we consider the first term inside the curly brackets. Here we are 
allowed to replace the spherical Bessel function with his leading term in the small argument 
expansion getting 



1 dse~ lAtS sm(sk F r) ~ - ——-L. e l 4A t . ( 2 1) 



3 ./ 2 2 V A t i 3 3 

Hence our statement follows. 

In Figs. H]-[6]we display the modulus of U°(x, 0) together with its real and imaginary parts 

for a few values of A t . We observe that for small A t a diffraction pattern emerges as in the 

case of the single-particle propagator: indeed, from the Heisenberg principle, here the energy 

is considerably spread out and, as a consequence, a wave packet can be set up which vanishes 

at fixed positions selected by the medium. As A t increases this pattern is washed out until 

for very large A t it turns into an almost uniform behaviour. From the formal point of view 



this evolution reflects the fact that for small time differences the particle-hole propagator 
is essentially imaginary, just as happened for the one- hole propagator. This imaginary part 
oscillates with the distance and hence the diffraction pattern follows. However, as the time 
difference increases a real part develops, also with an oscillatory behaviour, but with zeros 
displaced with respect to those of the imaginary part. Hence the zeros of the diffraction 
pattern are lifted up until at large At the modulus of U°(x, 0) becomes uniform in space. 

A further important feature of U°(x, 0) relates to its rapid decrease occurring in the range 
< kpr < 4: we have numerically checked that in this domain all the terms in Eq. ( T20|) 
contribute by roughly the same amount whereas for larger kpr the integrals with the variable 
running in the intervals (0, 1) and (1,2), where the Pauli correlations are operative, essen- 
tially drop out. We thus conclude that are these correlations which damp the propagation 
of a density disturbance in the system. 

We turn now to the evaluation of the remaining integrals with respect to the variable q 
to complete our task of obtaining an analytic expression for n°(x, 0). For this purpose we 
introduce the dimensionless quantities 

kpr 



2^/A t 

and the function 



and p = a/ A t (22) 



rb e iy 2 

g(a,b,c) = I dy— — . (23) 



The latter can in fact be expressed in terms of the Meijer-G functions [8j, although the 
resulting formulas are quite cumbersome and hence we prefer to use directly the definition 
expressed by Eq. f)23p . Even so the remaining g-integrations are, unfortunately, given by 
expressions which are far from simple; we report these in Appendix [C] for completeness. 

To pave the way to Appendix [0 here we simply rewrite Eq. fl20l) in terms of the variables 
in Eq. (j22|) . It becomes 

Il a (x,x') = -inl- — 3 / dse- ip2s \\n(2pzs)j 1 {2p 2 s) 

P J 2 

~ -inl- I dse~ ip2s2 s 2 for p -> (24) 
2 J 2 

U° b (x,x') = -ir%~ j\se-^ s2 sm(2pzs) 3l (2p 2 s(s-l))(s-l) 2 
49 

- ~ in o^ forp^O (25) 
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U°Jx,x') = -inl— 4 I dssm(2psz) \e~ ip2s2 cos(2ps(l - s)) - e - 2 '"'' v '' 
16z p 5 I 



x 



1 

2/7s 



cos(p 2 s(s - 2)) + 7rT - sin(p 2 s(s - 2)) 
21 1 



- m l^jA forp^O (20) 
9 1 fcos(4pz)-l 3 f 1 /o„„x-a<pi» a „• /o „2 2N 1 / 1 



n°(x,x') = -inl — i_/ cos ( 4 ^) ^ 3 + /* ^ S i n ( 2 p^)e- 2 ^ s sin(2p 2 s 2 )- 
loz p & [ 8zp d ,/ s 



2sp 



2 



,2 



9 z 

— ""n — 7 for p ^ (27) 

Concerning this all important issue, referred to as second kind scaling, much light on it 
is shed by the analysis of II (x, 0) in space-time coordinates. Here one realizes that just as 
G°(x, x') turned out to be proportional to n (see Eq. (Tjl|)), U°(x, 0) is proportional to Uq as 
expected. Then, when these density factors are divided out in both propagators one finds 
that G°(x,x') and n°(x, 0) scale at any A t providing that the space-time coordinates are in 
turn rescaled in terms of the Fermi momentum and frequency. This is clearly illustrated in 
Figs. HI and [H] where the spatial behaviour of the n°(a;, 0) associated with three different 
values of kp is displayed as a function of kpr. The curves are seen to coincide at A t = 0.8, 
1 and 1.2, in accord with Eq. (120]) . which transparently exhibits the fcp-scaling property. 

The above results concerning U°(x, 0) should, however, be viewed with some care. This 
is because the search for scaling in kp at vanishing At is impossible owing to the divergence 
discussed above which signals that the theory is ill-defined. Actually all of the terms in 
Eqs. ( I20l) diverge at A t = 0, with the exception of the second one. 



V. THE IMPACT OF A MORE REALISTIC MOMENTUM DISTRIBUTION 

Here we explore the impact of interactions among the constituent fermions on the pair 
correlation function and the CSR. This we do in a schematic frame that should, however, 
capture some of the relevant physics. We assume for the momentum distribution n(k) the 
expression 

n(k) = 6{k F - k){l - a-j) + 9(k - k F )p 1 e~ (32{ ^~ 1) . (28) 
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The four parameters (indeed also k F should be viewed as such) entering in Eq. ([28]) must 
satisfy the normalization condition 



7T Z 



i-? + ^(/?l + 2/? 2 + 2) 



(29) 



.3 5 ' PI 

where uq is the experimental constant density of the system. Notice that for 1 — a — fi\ the 
Fermi system becomes superconductive (the Fermi surface disappears) whereas for 1—a > Pi 
the discontinuity at the Fermi surface remains, as it should for a normal Fermi system 
according to the Luttinger H[ theorem. This implies that both a and Pi should be positive. 

We now choose as an example nuclear matter (in this case the left-hand side of Eq. (|29|) 
should be multiplied by 2 to account for the isospin degeneracy) where one has no = 0.17 
fm -3 . We then display in Fig. [7] the n(k) for nuclear matter for a specific choice for the four 
parameters, chosen to fulfill both Eq. (T2"9~|) and the above-mentioned constraint. The tail 
at large momenta is evident and one finds that the new Fermi momentum turns out to be 
kp = 1.54 fm -1 , namely larger that the one associated with the non- interacting case. 

The pair correlation function for the momentum distribution in Eq. (128]) is then easily 
computed and reads 

1 



i r 3 



2{k F r 



ji{k F r) - 



a 



(k F r) 4 



2>{{k F r) 2 - 2) sm(k F r) 



-k F r{{k F r) 2 — 6) cos(k F r)) + Pi 
Pl + P 



k F r 



(k F r) 2 + p\ 



(30) 



x ( sin(^) ( A + «±£ - l) + cos(^) (k F r + 2^±S) ) 



(k F r) 



k F r 



Using the same values of the parameters as for the momentum distribution shown in Fig. [71 
this pair correlation function is displayed in Fig. [8J where it is compared with that of the 
pure Fermi gas from Eq. (|11|) . What clearly appears in the figure is the marked difference 
between the two correlations functions at short distances, while they practically coincide at 
large distances: this behaviour nicely illustrates the role of the short-range correlations. 

Finally, we compute the CSR using Eq. (1301 . Although also in this case the calculation 
can be analytically performed using complex coordinates, as was previously done in the 
non-interacting situation, the resulting expression turns out to be very cumbersome; hence 
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we resort to the numerical evaluation of the formula 

S(q) =Z-v^J dfe-^g 2 (r) , (31) 

the function g(r) being defined in Eq. (|30|) . 

The outcome is shown in Fig. [9] where it is clearly seen that results from Eq. fl3T|) coincide 
with those from Eq. (fl3l) at large q (say q > 4 fm _1 ), as they should, since in this domain 
of momenta the associated wavelengths are so small that the system appears to the probe 
as a collection of uncorrelated fermions. On the other hand for, say, 2.5 < q < 4 fm -1 the 
two sum rules differ substantially due to the action of the correlations among the fermions. 
Finally, for smaller q this difference tends to disappear as both sum rules should go to 
zero when q vanishes. Noteworthy is that the sum rule arising from Eq. (I3T]) reaches the 
asymptotic value 1 from below. 

VI. CONCLUSIONS 

In this work we have deduced expressions for the one- and two-particle propagators 
as functions of space-time coordinates, focusing on the hole propagator and the density- 
density correlator. To our knowledge these expressions were not previously available in the 
literature. We find that both propagators have infinities at equal time, a problem that is 
being addressed in other work. Next we have explored the asymptotic space-time behaviour 
of both Green functions and have found a transition from a diffractive regime to a uniform 
one as the time difference between the fields (for the one-particle propagator) or between the 
densities (for the two-particle propagator) grows. From the formal point of view this relates 
to the fact that both propagators for zero time difference are purely imaginary, but then the 
real parts start to develop as the time difference grows. Concerning the dependence upon 
kp, our analysis shows that both G° and n scale once appropriate measures for space and 
time are chosen. This outcome goes in parallel with the situation in frequency-momentum 
space. However the divergence affecting II at At = prevents the analysis of second-kind 
scaling for this propagator at vanishing At. 

We have found that the key ingredient contained in the propagators is the momentum 
distribution n(k). For example the Coulomb sum rule can be directly expressed in terms 
of the pair correlation function, which can, of course, be obtained once n(k) is known. We 
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have explored the consequences arising from a modification around the Fermi surface of 
n(k) away from the pure Fermi gas result, finding that this induces striking effects in the 
pair correlation function at short distances. This in turn leads to major differences in the 
Coulomb sum rule for momenta between about 2.5 and 4 fm -1 for the model used in the 
present study, suggesting that it would be interesting to explore the responses of our infinite 
Fermi system to external probes, employing in the calculation of the polarization propagator 
our modified momentum distribution function given in Eq. (128]) . This, of course, is meant 
to account for the correlations (possibly of short range) among the fermions. 

A further issue, in some sense complementary to the above one, deserves consideration. It 
relates to the passage from an infinite to a finite system, the latter obviously of concern for 
nuclear physics. As a first approximation this transition can be accomplished by accounting 
for the modification of the density of the single-particle states induced by the presence of 
the system's surface according to the prescription given by Feshbach 12j. In a preliminary 
investigation we have done so and found, as expected, that the impact of the surface in 
coordinate space on the system's response is only felt at low momenta. More specifically, 
the confinement on the one hand entails oscillations of the system's density in coordinate 
space, on the other enlarges the momentum distribution n(k) to momenta greater than those 
of the corresponding infinite system with equal density. At the same time it digs a hole at 
low k in the momentum distribution. Thus confinement (in leading order) and short-range 
correlations appear to work in the same direction at large momenta. The disentangling of 
the interplay between the two effects is a problem we intend to address in forthcoming work. 
We will also explore in depth how the scaling of first kind is reflected in the space-time 
coordinates. 

Finally, since much of the physics we are addressing here occurs at large momenta, it 
is imperative to extend the present treatment to the relativistic context, with the goal of 
providing a comparison between the results one would obtain using the present approach 
and those already obtained in other approaches to electroweak superscaling. For instance, 
it will be of interest to answer the question: How will scaling in kp of the imaginary part 
of the polarization propagator occurring for all values of space-time coordinates in the non- 
interacting, homogeneous case be affected (and eventually disrupted) by confinement and 
short-range correlations in both the non-relativistic and relativistic context? 
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APPENDIX A 

In this appendix we analyze the behaviour of the propagator G°(x,x') for very small A t . 
For this purpose we use the well-known asymptotic expansion of the error function 

w (2m-l)! 



1 _ 2 

erf(z) = 1 - —=-e z 



(2z 2 ) m 

m=l v ' 



(Al) 



Then from the above, after some algebra, in leading order of A t one obtains for the third 
term on the right-hand side of the propagator in Eq. (jl]) 

1 kplx — x'\ ( e ik F \x-x'\ e ~ik F \x-x'\ 



2 2A t \k F \x - x'\ + 2A t k F \x-x'\-2A i 

1 . /, ,-, cos(k F \x - x'\) . . n . 

sm(k F \x-x\) p— i — — — (A2) 



2iA t k F \x — x 

which exactly cancels the second term of the right-hand side; hence Eq. (J2J) follows. 

APPENDIX B 

Here we compute the integral in Eq. (1121) working in the complex z-plane. Setting q 
q/kp, we have 

'sinz cosz\ 2 sm(qz) 



f°° 1 f°° (i 

J dzji{z)j {qz) = - J dz\^- 



z 2 z J qz 



1 f°° dz 

= — — (sin z — z cos z) 2 sin(gz) . (Bl) 

2<? J -OO Z 

The integrand, which behaves as ~ z 2 for z — > 0, is a regular analytic function; hence the 
integration path along the real axis can be deformed by inserting a very small semicircle 
going around the origin from below. This we indicate with the symbol -f. Closing the 
integration path along a large semicircle in the Im z > region, we thus get 

/ = 4f ^ { (1 - %zfe^ z - (1 - izfe^ z + (1 + iz) 2 e^- 2 > 

°Q J -oo Z ^ 

-(1 + izfe- l{ ^ 2)z - 2(1 + z 2 )(f qz + 2(1 + z 2 )e- mz } , (B2) 
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where the fourth and the sixth term on the right-hand side do not contribute because for 
these the integration path has to be closed in the Im z < domain where no singularities 
exist. It is then convenient to split Eq. (lB2p into two pieces according to 



I = h + h 



(B3) 



with 



and 



h 



dz 



1 - iz) 2 e i( 3 +2 > - 2(1 + z 2 )e 



iqz 



(B4) 



£^^(l+i2)2 e i(?-2)* if g > 2 



(B5) 



v / I - £^ J(i _ lz f e ^-- q )z if g < 2 . 
The straightforward, while somewhat tedious, computation of the residues then yields 

/> = --{-(V-i 2)'-^w-! —(77 -i 2f-9^ ? 1- . (B6) 



and 



8q 1 4! 



^/i(g-2) 4 + |(g-2) 3 + i(g-2) 2 } if g > 2 



(B7) 



^ = +^{^(^-2) 4 + |(g-2) 3 + ^(g-2) 2 } ifg<2 . 
As a result of the cancellation occurring between Eqs. (IB6I) and (IB7I) one thus finds 



(B8) 



and 



I = I 1 + I 2 = Q if 5 > 2 



I = I 1 + I a = |(i-* + i*») ifg<2 



(B9) 



(BIO) 



APPENDIX C 

In this appendix we perform the integrals in Eq. (1191) . We focus on the first piece, since the 
second one is immediately derived once the first is known. For these pieces three angular 
integrations are trivial whereas the fourth should be approached with care owing to the 
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presence of the ^-function, leading naturally to the splitting of the integral over the modulus 
of the vector q into three pieces: 



n°(x,o)0(t a 
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The remaining angular integration is trivial and the integral over the modulus of the vector 
k, while somewhat cumbersome, can be performed. Introducing for convenience a time 
variable with the dimensions of a length squared, r = t x /m, one arrives at the expression 
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Then, for the first term in Eq. (1C2j) . the q- integration between Ik p and oo yields: 
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for the second term in Eq. (1C2j) the g-integration between kp and 2k p yields 
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for the term embodying the j\ and 
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for the other terms. Finally for the third term in Eq. (IC2I) the g-integration between and 
A;^ yields: 
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It is worth noticing that the above formulas embody the physics of diffraction (through the 
familiar erf) and the attenuation of the propagator (through the Meijer g-functions). 
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FIG. 1: The equal-time hole propagator normalized to the density no as given in Eq. ([7]) versus 
kpr (r is the modulus of the relative distance). 




(c) 

FIG. 2: Real part, imaginary part and modulus of G°(x,x') normalized to the density no for 
A t = 0.5 . 
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(c) 

FIG. 4: Real part (a), imaginary part (b) and modulus (c) of n°(x,0) divided by tiq and plotted 
versus kpr for At = 0.8 and three values of kp: 1.2, 1.36 and 1.5 fm _1 . Note that the same curve 
is obtained for any value of kp. 
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(c) 

FIG. 6: As for Fig. but now with A t = 1.2 . 
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FIG. 7: The momentum distribution of an interacting Fermi system as given by formula Eq. (|28|) 
of the text. The following values for the parameters have been chosen fcp = 1.54 fm _1 , a = 0.2 , 
(3\ = 0.4 and 02 = 4 . The non interacting case (dotted line) corresponds to kp = 1.36 fm -1 . 




FIG. 8: The pair correlation function 1 — ^g 2 (kpr) for a free (dotted line) and for an interacting 
(continuos line) Fermi gas as given by our model (formula (|30p of the text). The values of 
parameters are the same as in Fig. [7] . Note that the two curves refer to different kp . 
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FIG. 9: The Coulomb sum rule (see Eq. (|3ip of the text) for a free Fermi gas (dotted line) and 
for a correlated one according to our model (continuous line). The parameters are the same as in 
Fig. [7] and Fig.®. 
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